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Dissociation of accumulated genetic risk and disease 
severity in patients with schizophrenia 

S Papiol 1 ' 2 ' 3 , D Malzahn 4 , A Kastner 1 , S Sperling 1 , M Begemann 1 , H Stefansson 5 , H Bickeboller 4 , K-A Nave 2 ' 3 and H Ehrenreich 1,3 

Genotype-phenotype correlations of common monogenic diseases revealed that the degree of deviation of mutant genes from 
wild-type structure and function often predicts disease onset and severity. In complex disorders such as schizophrenia, the 
overall genetic risk is still often >50% but genotype-phenotype relationships are unclear. Recent genome-wide association 
studies (GWAS) replicated a risk for several single-nucleotide polymorphisms (SNPs) regarding the endpoint diagnosis of 
schizophrenia. The biological relevance of these SNPs, however, for phenotypes or severity of schizophrenia has remained 
obscure. We hypothesized that the GWAS 'top-10' should as single markers, but even more so upon their accumulation, display 
associations with lead features of schizophrenia, namely positive and negative symptoms, cognitive deficits and neurological 
signs (including catatonia), and/or with age of onset of the disease prodrome as developmental readout and predictor of disease 
severity. For testing this hypothesis, we took an approach complementary to GWAS, and performed a phenotype-based genetic 
association study (PGAS). We utilized the to our knowledge worldwide largest phenotypical database of schizophrenic patients 
(/?>1000), the GRAS (Gottingen Research Association for Schizophrenia) Data Collection. We found that the 'top-10' GWAS- 
identified risk SNPs neither as single markers nor when explored in the sense of a cumulative genetic risk, have any predictive 
value for disease onset or severity in the schizophrenic patients, as demonstrated across all core symptoms. We conclude that 
GWAS does not extract disease genes of general significance in schizophrenia, but may yield, on a hypothesis-free basis, 
candidate genes relevant for defining disease subgroups. 
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Introduction 

In complex disorders lacking clearly identifiable disease- 
causing factors, such as schizophrenia, the overall genetic 
risk is often >50% 1 but genotype-phenotype relationships 
are obscure. Recent genome-wide association studies 
(GWAS) on schizophrenia, building on very large cohorts of 
cases and controls, have uncovered and replicated a risk for 
several single-nucleotide polymorphisms (SNPs). Among the 
10 'top hits', reaching genome-wide significance in different 
populations after multiple-testing correction (even though with 
low odds ratios), are markers in the major histocompatibility 
complex region and in TCF4, ZNF804A and NRGN genes. 2-6 
The biological relevance of these SNPs, however, for 
phenotypes or severity of schizophrenia has remained 
unclear. 

All current GWAS data on schizophrenia rely on endpoint 
diagnosis only and do not allow for genotype-phenotype 
correlations. In an approach complementary to GWAS, we 
employed the 'top-10' schizophrenia-associated SNPs that 
have emerged as genome-wide significant from several 
GWAS, and explored their potential contribution to the 
disease phenotype, including positive and negative symp- 



toms, cognitive deficits and neurological signs (including 
catatonia), and age of onset of the prodrome as develop- 
mental readout and predictor of disease severity. For this 
phenotype-based genetic association study (PGAS), we 
utilized the GRAS (Gottingen Research Association for 
Schizophrenia) Data Collection. 7 ' 8 

Materials and methods 

Schizophrenic patients. The GRAS study was approved 
by the Ethics Committees of the Georg-August-University of 
Gottingen and of participating centers, and comprises at 
present 1041 patients with confirmed Diagnostic and Statis- 
tical Manual of Mental Disorders, fourth edition (DSM-IV) 9 
diagnosis of schizophrenia (82.2%) or schizoaffective 
disorder (17.8%), examined between 2005 and 2010 in 23 
collaborating centers all over Germany (Supplementary 
Table 1). 7 ' 8 

Healthy subjects. Healthy subjects for the case-control 
study were blood donors (n=1144), recruited according to 
national guidelines for blood donation. 7 Comparable to the 
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patient population (Caucasian 95.5%; other ethnicities 1.8%; 
unknown 2.7%), almost all controls were of European- 
Caucasian descent (Caucasian 97.8%; other ethnicities 2%; 
unknown 0.2%). 

Phenotyping. Comprehensive interviews, testing and 
clinical ratings were conducted by an invariable team of 
trained examiners (psychiatrists and psychologists) using the 
'GRAS Manual'. 7 ' 8 Diagnoses of schizophrenia/schizo- 
affective disorders were based on the structured clinical 
interview for DSM-IV (SCID), 10 substantiated by information 
from medical records/discharge letters of all the patients. 
Psychopathological state, symptom severity and functional 
outcome were evaluated by clinical ratings (positive and 
negative syndrome scale (PANSS) 11 and global assessment 
of functioning (GAF). 9 ) Neuropsychological testing (subtest 
3 of Leistungsprufsystem, 12 Trail-Making Test, 13,14 Verbal 
Learning and Memory test 15 ) including pre-morbid intelligence 
(Mehrfachwahl-Wortschatz-lntelligenz test-B 16 ) and neuro- 
logical examination (Cambridge Neurological Inventory 17 ) 
yielding respective composite scores, are described in detail 
elsewhere. 7 ' 8 

Genotyping. SNP genotyping was performed with 
SimpleProbes (TIB Molbiol, Berlin, Germany) on 
LightCycler480 (Roche, Mannheim, Germany). All markers 
fulfilled the Hardy-Weinberg equilibrium. 

Statistical methods. PGAS: Phenotype data were 
standardized to be normally distributed with expectation 
zero and variance one, and presented such that higher 
values always indicate better performance. Expected values 
of univariate and multivariate phenotypes were analyzed with 
linear models (Rv2.12.0), covariate adjusted and tested for 
additive effects of GWAS risk alleles with single-locus 
models (Table 1) and with a 10-loci model (Table 2). The 
latter simultaneously estimates for all the 10 loci regression 
coefficients with number of respective GWAS risk alleles. 
To assess statistical significance, the commonly used 
Bonferroni correction was employed; however, raw 
P-values are given. 

Results 

First, we proved in a case-control study (cases n=1041; 
healthy controls n=1144) that the GRAS population (for 
patient characteristics see Supplementary Table 1 ) provides a 
'genetic data matrix' that essentially replicates the GWAS 
results (Table 1, upper part). In fact, screening of the 10 
genome-wide hits resulted in a significant association of 
markers rs6913660, rs1 321 1507 and rs3131296 in the major 
histocompatibility complex region (chromosome 6), and 
marker rs2312147 in chromosome 2 (VRK2 gene) with an 
increased risk for the disease. Due to the smaller sample size 
compared with the GWAS studies, leading to lower power 
(power 70-90% for the SNPs found to reach allelic P-values 
< 0.05, average power of 43% over all markers), not all SNPs 
turned out to be significantly associated with the schizo- 
phrenia risk. Nevertheless, all of them exhibit the same 
direction of association as reported in GWAS. 2-6 In fact, the 



GRAS sample has been included in a recent large GWAS 
follow-up study. 6 

We next conducted PGAS single-locus and PGAS multi- 
locus quantitative association analyses. In both the proce- 
dures we searched for the hypothesized association of 
markers with the lead symptoms of schizophrenia including 
developmental readouts/disease severity, both in the form of 
a composite construct (multivariate phenotype) and of its 
components separately (individual phenotypes), that is, 
positive and negative symptoms (PANSS), cognitive deficits 
(composite of executive function, reasoning, verbal learning 
and memory) and neurological signs (including catatonia), 
and/or with the age of onset of the disease prodrome. 
We further hypothesized that, if the markers were disease 
specific, they should not be associated with the schizo- 
phrenia-unrelated, general disease control variables. 
Hence, such variables referring to symptoms that are not in 
any way specific for schizophrenia were also included 
PANSS general psychopathology (depression, anxiety and 
others), global assessment of functioning and pre- 
morbid intelligence (basic cognitive capabilities of an indivi- 
dual before disease onset) (Tables 1 and 2; Supplementary 
Table 2). 

With >1000 patients, the single-locus analysis found 
none of the schizophrenia phenotypes, neither individual 
nor multivariate, significantly influenced by any of the markers 
after multiple-testing correction by Bonferroni (Table 1). 
We thus wondered whether — instead of single genotypes — 
the 'genetic load' of a patient, that is, the accumulation of risk 
genotypes, would give us a clearer signal regarding the 
contribution of genetic risk to disease phenotype and severity. 
As illustrated in Figure 1 , the composite severity score, built on 
the basis of the five core symptom variables (see inset), 
shows an essentially equal distribution of schizophrenia 
severity from the lowest to the highest genetic load group. 
As more symptoms do not necessarily reflect more 
severe disease and certain symptom groups may better 
associate with genetic load, all five variables were considered 
separately, too. But again, for the 10 top hits of GWAS, none 
of the schizophrenia symptom groups was dependent 
regarding its severity on the increase in the 'genetic load' 
(Supplementary Table 2, and Supplementary Figure 1). This 
becomes evident when the composite score data are 
presented alternatively as percentage severity over all the 
genetic load groups; there is no tendency of an increase 
in severe cases with an increase in the number of risk 
markers (Supplementary Figure 2). Interestingly, grouping the 
distribution of accumulated risk genotypes in the Icelandic 
GWAS sample (n = 582 schizophrenic individuals) yields a 
pattern of overall risk distribution that is comparable to that 
of the GRAS database (Supplementary Figure 3), further 
supporting the validity of our sample for the PGAS approach to 
the GWAS hits. 

We now pursued the question whether better insight into the 
phenotypical contribution of the GWAS hits could be obtained 
by performing a multilocus additive joint model analysis. This 
analysis weighs each marker according to its estimated 
relative importance as a risk genotype. Also this approach 
failed to uncover statistically significant associations after 
multiple-testing correction. Close to significance are two 
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Table 1 GWAS replication and PGAS single-locus quantitative association analysis 



Region/ neighboring gene 
(Chromosome) 


MHC 
HIST1H2BJ 
(Chr.6) 


MHC 
PGBD1 
(Chr.6) 


MHC 
NOTCH4 
(Chr.6) 


NRGN 
(Chr.11) 


ZNF804A 
(Chr.2) 


MHC 
PRSS16 
(Chr.6) 


VRK2 
(Chr.2) 


TAOK2 
(Chr.16) 


TCF4 
(Chr.18) 


TCF4 
(Chr.18) 


SNPID 


rs69 13660 


rs1 321 1507 


rs3131296 


rs1 2807809 


rs1 344706 


rs6932590 


rs2312147 


rs4583255 


rs9960767 


rs4309482 






C 
C 


73.4/65.8 


T 
T 


86.2/81.9 


G 
G 


79.3/74.4 


T 
T 


69.7/66.7 


T 
T 


38.2/35.2 


T 
T 


61.0/57.1 


C 
C 


40.8/35.3 


T 
T 


33.8/31.5 


A 
A 


88.3/89.0 


A 
A 


37.7/35.6 


Genotype 
Case(%) / Control (%) 


c 

A 


24.3/31.3 


T 

C 


13.2/17.5 


G 
A 


19.6/24.1 


T 

C 


27.8/30.1 


T 

G 


49.0/49.3 


T 

C 


33.7/37.4 


C 
T 


43.3/49.4 


T 

C 


47.7/47.5 


A 
C 


11.4/11.0 


A 
G 


46.3/47.2 






A 
A 


2.3/2.9 


C 

c 


0.6/0.6 


A 
A 


1.1/1.5 


C 

c 


2.6/3.2 


G 
G 


12.9/15.4 


C 

c 


5.3/5.5 


T 
T 


15.9/15.3 


C 

c 


18.5/20.9 


C 
C 


0.3/0.1 


G 
G 


16.0/17.3 


P value 


0.0007 




0.0195 


0.0255 




0.3036 




0.1562 


0.1707 




0.0130 


0.3054 




0.5058 




0.5301 


Allelic 


C 


85.5/81.4 


T 


92.8/90.6 


G 


89.1/86.5 


T 


83.5/81.8 


T 


62.6/59.9 


T 


77.8/75.8 


C 


62.5/60.0 


T 


57.6/55.3 


A 


94.0/94.4 


A 


60.9/59.2 


Case(%) / Control(%) 


A 


14.5/18.6 


c 


7.2/9.4 


A 


10.9/13.5 


c 


16.5/18.2 


G 


37.4/40.1 


C 


22.2/24.2 


T 


37.5/40.0 


C 


42.4/44.7 


C 


6.0/5.6 


G 


39.1/40.8 


P value 


0.0004 


0.0094 


0.0081 




0.1264 




0.0695 


0.1122 


0.1042 






0.5568 




0.2628 


OR (Risk Allele) 
(95% CI) 


(1 


1.34(C) 
.14, 1.59) 


(1 


1.34 (T) 
07, 1.68) 


1.28 (G) 
(1.06, 1.55) 


1.13 (T) 

/n OR A QO\ 

(u.yo, \ .oo) 


(0 


1.12 (T) 
.99, 1 .27) 


1.12 (T) 
(0.97, 1.30) 


(0 


1.11 (C) 
.98, 1.25) 


1.10 (T) 
(0.97, 1.24) 


1.08 (C) 

{U.OO, 1 .4 I ) 


1.07 (A) 
(0.95, 1.21) 


OR (Risk Allele) 

(95% CI) 
in previous GWAS 


1.15(C) 
(1.10, 1.21) 


1.24 (T) 
(1.16, 1.32) 


1.19(G) 
(1.13, 1.25) 


1.15 (T) 
(1.10, 1.20) 


1.10 (T) 
(1.07, 1.14) 


1.16 (T) 
(1.11, 1.21) 


1.10(C) 
(1.06, 1.13) 


1.08 (T) 
(1.05, 1.11) 


1 .23 (C) 
(1.15, 1.32) 


1.09 (A) 
(1.06, 1.12) 


Multivariate 


Allele Effect 


0.0054 


0.0286 


-0.0697 


0.0979 




0.0174 


0.0027 


0.0263 


0.0064 


0.1158 




0.0066 


phenotype 3 


P value 


0.8837 


0.5731 


0.0987 


0.0057 


0.5280 


0.9322 


0.3107 


0.8072 


0.0346 


0.8014 


PANSS 


Allele Effect 




0.0014 


0.0736 


-0.0245 


0.0317 




0.0043 


0.0157 


0.0278 


0.0522 


0.1054 




0.0380 


positive 


P value 


0.9811 


0.3767 


0.7237 


0.5838 


0.9233 


0.7597 


0.5121 


0.2228 


0.2396 


0.3779 


PANSS 


Allele Effect 


0.0004 


0.1260 


0.0542 


0.1572 


0.0066 


0.0025 




0.0031 


0.0813 


0.1673 


0.0587 


negative 3 


P value 


0.9951 


0.1398 


0.4471 


0.0081 


0.8862 


0.9624 


0.9441 


0.0653 


0.0677 


0.1870 


Cognitive 


Allele Effect 




0.0521 




0.1136 


0.0441 


0.1301 


0.0003 




0.0133 




0.0342 


0.0433 


0.0851 




0.0298 


score 3 


P value 


0.2802 


0.0857 


0.4225 


0.0048 


0.9932 


0.7459 


0.3169 


0.2063 


0.2382 


0.3870 


Total CNI 3 


Allele Effect 


0.0458 




0.0135 


-0.1239 


0.1626 




0.0213 


0.0229 




0.0137 


0.0410 


0.0436 




0.0030 


P value 


0.4307 


0.8655 


0.0639 


0.0033 


0.6197 


0.6401 


0.7373 


0.3190 


0.6070 


0.9425 


Prodromal 


Allele Effect 


0.0148 


0.0462 




0.1974 


0.0530 




0.0453 


-0.0204 


0.0941 


-0.0976 


0.1480 


0.0097 


Onset 


P value 


0.8256 


0.6259 


0.0103 


0.4174 


0.3683 


0.7179 


0.0428 


0.0387 


0.1472 


0.8378 


PANSS 


Allele Effect 




0.0025 


0.0656 


-0.0040 


0.1392 




0.0350 




0.0026 


0.0340 


0.0678 


0.1334 


0.0244 


general 


P value 


0.9679 


0.4506 


0.9562 


0.0213 


0.4586 


0.9604 


0.4433 


0.1290 


0.1524 


0.5885 


GAF 3 


Allele Effect 


0.0035 




0.0382 




0.1233 


0.1147 




0.0139 


0.0378 


0.0157 


0.0307 


0.1199 


0.0196 


P value 


0.9552 


0.6569 


0.0858 


0.0558 


0.7658 


0.4743 


0.7204 


0.4891 


0.1960 


0.6616 


Premorbid 


Allele Effect 




0.0435 




0.0434 


0.0309 


0.0542 


0.0242 




0.0020 


0.0612 




0.0824 


0.0100 




0.0067 


Intelligence 13 


P value 


0.4770 


0.6083 


0.6622 


0.3589 


0.6008 


0.9700 


0.1617 


0.0589 


0.9134 


0.8786 



Abbreviations: CNI, Cambridge neurological inventory; CI, confidence interval; Chr, chromosome; GAF, global assessment of functioning; GWAS, genome-wide 
association studies; MHC, major histocompatibility complex; OR, odds ratio; PGAS, phenotype-based genetic association study; PANSS, positive and negative 
syndrome scale; SNP, single-nucleotide polymorphism; SNP ID, SNP identifier. 

GRAS sample of schizophrenic patients, n=1041; healthy control sample, n=1144. Upper part shows the case-control genetic association study essentially 
replicating previous GWAS results. SNPs are presented from left to right in the order of OR (odds ratio). Pearson / 2 -test and Fisher's exact test (both two sided) were 
used for genotypic and allelic comparisons, respectively. Lower part (PGAS) gives additive effect per copy of GWAS risk allele on expected value of schizophrenia- 
relevant quantitative phenotypes. The multivariate phenotype combines five schizophrenia core features (positive PANSS, negative PANSS, cognitive score, total 
CNI and prodromal onset). All phenotypes were standardized to zero mean and variance one and presented such that larger values correspond to better performance 
(for this purpose, PANSS scores and total CNI were multiplied by -1). General PANSS, GAF and premorbid intelligence were included in the analyses as disease 
control variables. The estimate of allele effect is negative if carriers of a GWAS risk variant perform worse. It is positive, if carriers of a GWAS risk variant perform better 
with respect to the expected trait value for the schizophrenic sample. Allele effect size on mean trait is quantified relative to trait variability (standard deviation). P- 
values below 0.05 were highlighted for optical guidance but are not significant after multiple-testing adjustment. 
a Corrected for age: PANSS negative, cognitive score and total CNI (for separate analyses and within multivariate phenotype), GAF. 
b Corrected for language problems: pre-morbid intelligence (898 with no language problems and 108 with correction for language problems). 
Exploratory exclusion of non-Caucasian subjects from the GRAS sample (n = 48; 4.5%) did not qualitatively alter any of the main findings in this Table. 



associations, one between the marker rs3131296 in NOTCH4 
(chromosome 6) and prodromal onset (developmental 
readout), the other between rs1 2807809 in chromosome 11 
(near NRGN gene) and severity of the neurological signs 
(Table 2). In fact, when ignoring the multiple-testing issue, the 
latter marker which has a relatively low odds ratio in our case- 
control study (OR 1.13 here and 1.15 in the GWAS study of 
Stefansson et al. 5 ), shows several association 'signals' that 
make it attractive for follow-up in future subgroup analyses. 
However, this marker did not associate with the cognition of 
schizophrenic subjects in a recent study. 18 

Discussion 

How can the overall negative result regarding the phenotypi- 
cal significance of basically all GWAS 'top-10' genotypes be 
explained? First, in clear contrast to monogenic diseases, 19 
the genetic risk for schizophrenia may not simply be reflected 



by phenotypical disease severity or by core symptoms of the 
disease. Here, an array of environmental risk factors that 
cannot easily be controlled for might also have a modulating 
role. 20 Second, the genotype -to phenotype translation may 
only be visible and valid for a relatively small subgroup of 
individuals, but still leads to significant (even though low) 
genetic risk odds ratios in very large GWAS samples. 
Third, one risk genotype may partly 'neutralize' another one, 
resulting in risk reduction upon combination rather than 
accumulation of the genetic load. This latter, seemingly 
paradox interaction is supported by the observation that the 
few effects on phenotypes found here in marker rs1 2807809 
(chromosome 11; near NRGN) unexpectedly go into the 
opposite direction (risk genotype shows less severity). Fourth 
and finally, we cannot rule out that for some analyses the 
GRAS sample may not have enough power to detect the 
(certainly weak if any) phenotypical consequences of the 'top 
10' GWAS hits. In this context, however, the general question 



Translational Psychiatry 



GWAS SNPs and schizophrenia severity 

S Papiol etal 



4 



CO 

< 
o 

Q_ 



I 5 ? 



2p 



ft c 

So 



T- 

CO _ 
CO 3 s 



N 



5£ 



I. 

I— CO 1 

O v: 



a. 



SB* 
5o 



0) 

c 

0) 

o 

I! 

c c 
.2 o 

CD -C 



1 



CM 

T- 

CO 
CM 



CO 



1 

55 



I 
I 

CM 



CM 
CO 



CO 
CO 

T- 

o> 

CO 



cd 

o lo 
d d 

I 



CO <tf 

co m 

CM CM 
i- o 



00 CD 

co 

O 00 
O 00 

d d 

I 



i- 00 

'sf CD 
CM 'sf 
O CO 

d d 



^J- CM 

CM 00 

O LO 

O CD 

d d 

I 



CD CD 
O CD 
C\J 

o 

d d 

I 



cd I s - 
cd 

CD CM 
O i- 



^J- CD 
C\J CD 
O CD 



<D 0 

CD — 



CD 
CL 

>^ 
o 



CL 
0 

CO 



00 00 

C\J CD 

^ C\J 

O CO 

d d 

I 



o 

CM CD 



i- CM 

d d 



CD 'sf 

^- oo 

O CM 



CM ^J- 
CD CM 
O 
O CO 



CM O 
CD CO 



00 O 
CD CO 
CO LO 
O CM 

d d 

I 



CO CM 

oo in 

O CM 



o o o o 



i- 00 

oo 

CD CM 
O i- 



co h*. 
o 

O CM 



i- CO 
CD CD 
O CM 
O 00 

d d 

I 



00 LO 
CD CD 
ID 
O LO 

d d 



K o 
^- CD 
CM 

O ID 



CD CO 

o o 

O CD 
O CD 

d d 

I 



LO CD 
i- CD 

cd 
i- o 



o oo 

LO i- 

o 

d d 



^j- CD 
oo oo 
i- o 

i- CM 



o o o o 



i- LO 

o 

d d 



CD CD 
CD co 

^ CM 

d d 



oo oo oo oo 



LO CM 

CD O 

CM CD 

o ^j- 

d d 



^ 00 

o 

i- LO 

O 00 

d d 



CD CM 
CM 

O LO 
O CD 



CO CO 

cm 

O LO 
O CD 

d d 



CO CD 
LO ^J" 
CM ^J- 
O K 

d d 

I 



o o 

I 



^ CD 
CM CM 

o ^j- 
d d 

I 



CD O 
^ CD 
^ LO 

o ^t" 

d d 



CM 00 
i- CO 
O h~ 
O CD 

d d 

I 



i- CM 

*sf CD 

i- CM 

o 

d d 

I 



i- CM 

^J- LO 

O LO 

O CD 

d d 



CD O 
O 1^ 

co I s - 
o ^ 
d d 

I 



CD CM 
O) i- 
O CM 
i- o 

d d 

I 



LO LO 

o 

o i- 

d d 



CD co 

CO 00 
LO i- 

O LO 

d d 

I 



CM CD 
LO CD 
O CM 

d d 

I 



CO CM 
00 CD 
LO CD 
O i- 

d d do 



CD CD 
CM CD 
CM O 
O CD 



00 CD 
CD CO 
CO O 
O 

d d 



CD 'sf 

CO ^j- 

CM CD 

o 

d d 

I 



i- o 
CD i- 
00 i- 
O ^1" 

d d 

I 



CM 1^ 
CD CO 
i- CD 
O CD 

d d 



LO LO 
00 CM 
CD CO 
O CO 

d d 



LO CD 

i- co 
o 

d d 

I 



O CO 

i- CD 

I s - Lo 

o co 

d d 

I 



^ LO 
CM i- 
CM CD 
i- CM 

d d 



LO t- 

1^ 1- 

^ CM 



O CD O CD 

do do 



<D CD 
CD — 



^J" LO 
O LO 
O CD 
O CD 

d d 



CD LO 
CD 00 
CD CM 



CM CM 
CD 
O CD 

d d 

I 



0 

O 0 

ll 



CM 00 
1^ O 
CD i- 
O t- 

d d 



cm 

CM 
^J- CD 
-r- O 

d d 

I 



00 o 
CD CD 
CO 
O LO 

d d 

I 



O 0 

ll 



LO CM 
CO CD 
LO CO 
t- O 



i- o 

CO CM 
CO o 
CM O 



O O O O 



^J" CM 
00 O 

O CD 

do do 



O LO 
O CM 
CD i- 



LO 
i- CM 
CD O 
O LO 



co 

00 CO 
LO K 

O LO 



o o o o 



CD 0 

0 — 



(D 0 

0 — 



LO CM 
CM CD 
^- CD 
O LO 

d d 

I 



CD CD 
i- CD 
CD i- 
O 

d d 



CD CD 

o 

O CD 
O CD 

d d 

I 



0 

O 0 

ll 



00 i- 

^- CO 

LO LO 

i- o 

d d 

I 



CD CO 

cm 

i- o 

O CD 



O LO 
CD O 
CM CD 

o 

d d 

I 



0 

O 0 

ll 



CO 
CO 



< 

CL 



CO 
CD 
0 

c 
CO 
CO 

z 
< 

CL 



CD 

o 

o 



o 



0 

c 

0 
CD 

CO 
CO 

z 
< 

CL 



CD CD 
CO CD 

o 

O 00 

d d 

I 



CD 1^ 
O i- 
O CD 



CD CD 

co ^i- 

K CD 

o o 
d d 

I 



CD O 
K CM 
LO CD 

O T- 

d d 



CD 00 

^J- CD 
O LO 

d d 



LO i- 

CO 'sf 

co 
o 




CD CD 
CD CD 
^ CM 
O LO 

d d 



00 00 

^J- CM 

i- CD 

O 00 

d d 

I 



^J- CD 
^J- CO 

o 

d d 

I 



0 

0) 0 

ll 



< CL < CL <CL < Cl < Cl < Cl < Cl < Cl < Cl 



0 

o 
c 

0 
CD 

0 
c 

"D 

!q 
o 

« E 

LL . 

< 0 
(D a. 



CO 
CO 



■Q CO CD 0 
CO CD ^ C 
CD ^ 5 CD 

to o - O 

C _ CO 

It 

S B § # 

cO cO — -5 

E o O E 

O CO — CD 

S5f I 

'to 5 £ Z 
5 c co CO 



^lo « £ 



o 5 

w CO 
CD 3 

co 

SI 

t co 

„ CD 

>^ 

'■E 

CO CD 

° ^ 



cb DC 

§ 2 



CO 

< 
o 

Q_ 



c c c ^ — 

CD O) CO >, 

O O CO -c ~ 

5 § co § _J .g 

E o ^ CO CO CO 

5 "5 ? CO 

- 0 CD < 



£ £ 



E r 



?< ?§ 

) - . 9- 

: 8^ 3 

: — CO Q. 
- cd 2 co 

£ < s 

M- Q_ 

° CD O 

w > ^ 

^ CD 

) "-5 P. V. 



<; co 

Is 



~ cr 

CO CO 



co ^ 



o 

X 



- CD o CD 

) CD ° -° 

: t (C O 

? 'C ^ 



.g co 

.2 ^ 
o CD 



°r o a)- 

_ Q_ Q. — 

o> o co 

= N CD 



0 CD 

^ CD 
CO N 

B CD 
CO 5z CD 

E cd E 

111= 



CO - 
P CO 

4*1 

CH © - 
D) CO 

^ co e 

CD =5 
Q. CD CD 

S £ 



cd .2 Z 

151 



8 S| 

w n E 
$Z o 

10 Q CD 

-° r 8: 



: CO ° 
^ CO cp 



! 0 > 



CD ■ 



CD < 

_C - 

I — CD D) 

o E co 

_o CO CD > " 

"id W -r^ ^2 



co % 
Z a) . 

CL CL. 
O E 

f c ■> 



S € & 

!§« 
I si 

CO CL t= 

s ° = 

CO CD co 

E -c -o 

CL CD 

'1 (DC 



o © c 



"O c o 



C .Cl 



> CD S 

£ Q. E 

S ° ^ 

R N CD 

u r ^ 

CD O CO 



co o 



t^-o w 5 o 



r- CD 

o -o w 

v - '-<x CO 

§11 

£r c ^ 
o <b~ 

is? 

oZ c 
O)C0 OT ~ 

111 



O CD 



E 

CO 

O 



-> o 

: ^o^ 
> S 50 

i C 0) 

11 1 

< cO cO 
; 0 w 
) > CO 



0) r b 

ill 

-O ■ 0 

I s*" 

•° a. co 

^. CO > 
0 73 "O 

£ 2 § 

CO C 

0 "0 CD 

1 o E 
H E o 

r- C 0 
■-ON 

52 'co o 

CO CO 11 
>^ 2 ^ 
0 CD N 

2 2 =g 

CO 0 CO 
CLQ.-D 
0 ~ c 
co ^ as 
co c co 
^ I 2 

^ c co 

£ CO CD 
CD 



co i= co 

0 "C t 

£ B E 

. g 0 R -o 



CO 

co "O 

co 0 ^ 

CO O LO 

!o 

I If 

o ° °- 

S c E 
HT-c: co 
z -1 CO 

O ^ CO 

_ CD < 



o & 6S o £ 



0 

c £ 

v» CD p 

2^g 
Q q5 4= 



8 

© 0 a 

5 CL O 

0 52 »- 
o 2 ° 
c _ M_ 
0 £ "o 
C2)| a) 

||| 



CD _ i_i 

> ^ 0 

:-e -q 3 

c O 3 

CD c CO 

° V c 

° 2 i 

CD" CL CO 



P O CD 
0 0,^ 



CO 

00 p 5 

CD E CO 
0 © O 

C= -Q ' 

co 2 o 

CO CL c 
Z 0 *r 
< CD ° 
0_ CO C 
=3 O 

0 ?'co 



03 "5 ^ 

LL CO o 

< > ns * 



■c Q-^ 
0 0 cr 

"D > © 



ico ( 



CO ^2 ^2 



■2 .2 ^ CO 



o co < 



< > ? ( 



Translational Psychiatry 



GWAS SNPs and schizophrenia severity 

S Papiol et al 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 



Genetic load + 



Figure 1 Distribution of phenotype severity and cumulative genetic load with 
respect to the number of GWAS-identified 'top-10' risk SNP alleles in the GRAS 
population (bar graph). Phenotype severity is based on a composite score of the five 
core features of schizophrenia, displayed as inset (including intercorrelations 
between items). Score range in the GRAS sample is divided into three equal parts 
and ranked as mild, medium and severe disease phenotype. The blue line denotes 
the comparison of the risk SNP allele distribution in the healthy control sample. 

arises of how much clinical significance a genetic association 
may have if several thousand patients are needed to reveal a 
tiny effect on disease severity or phenotype. Thus, building on 
>1000 patients, one would at least have expected some 
more signals to pop up (that is, more nominal P-values around 
0.05 in the PGAS part of Table 1 ). In fact, the estimated power 
assessed according to Lettre ef a/., 21 who performed simula- 
tions in a similar context (sample size, normally distributed 
phenotype), amounts to overall around 80% at the a = 0.05 
level for the PGAS approach. The presented additive model, 
chosen to match our hypothesis of cumulated risk, performed 
similarly to a co-dominant one but appeared to be slightly 
more powerful. 

After all, we note that entirely different genotypes that were 
never found to be significantly associated with any schizo- 
phrenia risk in GWAS, still may profoundly modulate the 
schizophrenic phenotype, for example, of genes encoding 
neuregulin-1 , complexin 2 or COMT. 7 ' 22 ' 23 On the other hand, 
GWAS finds may be of general rather than disease-specific 
significance. Several studies have for instance suggested 
that schizophrenia and affective disorders are on a continuum 
of liability. Genetic linkage and association studies have 
proposed common disease loci for both the disorders 224 
Likewise, family studies show that first-degree relatives of 
bipolar patients have a higher risk for schizophrenia compared 
with first-degree relatives of healthy controls. 25 ' 26 Also other 
psychiatric diseases like alcoholism or major depression have 
been found to be associated with certain schizophrenia risk 
genes, for example, DISC1. 27 Thus, exploration of many other 
phenotypes available in the GRAS database, including 
candidate intermediate phenotypes 28 ' 29 or those reflecting a 
more dimensional approach to the disease, 30 might potentially 
be interesting. Purpose of the present study, however, was to 
cover the lead symptoms of schizophrenia in the first place. 

To conclude, GWAS approaches in diseases as complex as 
schizophrenia do not lead to the reconstruction of a 'common 



disease mechanism' or to the discovery of 'classical disease 
genes', as such genes obviously do not exist. What makes our 
study important for the clinician is that we can show, for the 
first time, that the combination of a whole battery of genetic 
pre-disposing factors (the 'top 10' GWAS finds in schizo- 
phrenia) in individual patients will not make their schizo- 
phrenic phenotype any different or worse than that of those 
patients who do not carry these genetic factors. Importantly, 
however, GWAS results may guide, on a hypothesis-free 
basis, to the identification of totally unexpected candidate 
genes involved in certain disease aspects in subgroups of 
patients, as they can be defined by PGAS. In order to get 
closer to understanding the disorders as complex and 
heterogeneous as schizophrenia, GWAS and PGAS will have 
to go hand in hand. 
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